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We report scanning Hall probe microscopy ol highly underdoped superconducting YBa2Cu306+^ 
with T c ranging Irom 5 to 15 K which showed distinct flux bundles with less than one superconducting 
flux quantum ($0) through the sample surface. The sub-"I>o features occurred more frequently for 
lower T c , were more mobile than conventional vortices, and occurred more readily when the sample 
was cooled with an in-plane field component. We show that these features are consistent with kinked 
stacks of pancake vortices. 



PACS numbers: 74.25.Ha, 74.25.Qt, 74.72.Bk 

I. INTRODUCTION 

Cuprate superconductivity occurs mainly in the ab di- 
rection on the CuC>2 planes. This quasi-two-dimensional 
nature manifests itself in the anisotropy between the o 
axis penetration depth (A c ) and the in-plane penetration 
depth (X a b)- Clem 1 showed that in highly anisotropic 
layered superconductors a c-axis vortex can be viewed as 
a stack of magnetically coupled "pancake" vortices, one 
in each layer. This formulation suggested the possibil- 
ity for novel vortex behavior and has become a major 
part of the phenomenological understanding of cuprate 
superconductors. However, direct observations of sepa- 
rated pancakes or pancake stacks have been rare^ 

In this paper, we present scanned probe microscopy 
of magnetic flux in highly underdoped YBa2Cu306-f:r 
(YBCO) single crystals. We observed nearly isolated flux 
features with less than one flux quantum ($0 = hc/2e = 
20.7 G fim 2 ) through the sample surface, which we call 
"partial vortices" . A model of separated pancake vor- 
tex stacks, similar to a kinked structure suggested by 
Benkraouda and Clem^ but with a more important role 
for pinning, agrees well with our observations. 

Non-quantized flux in superconductors has been ob- 
served experimentally, arising for different reasons. Geim 
et al£ observed non-quantized flux penetration in meso- 
scopic thin film samples of aluminum due to two effects: 
the proximity of the vortices to the sample edge, and a 
surface barrier to flux penetration. Sub-$o flux has been 
imaged in YBCO thin films along grain boundaries sepa- 
rating regions of the crystal rotated 45° about the c-axis 
due to the d-wave symmetry of the pairing state in com- 
bination with facets along the grain boundary*^ For our 
measurements discussed here, the vortices were far from 
the edge and there were no rotations of the crystal axes 
aside from 90° twinning, so neither of these mechanisms 
applies. Our observations of seemingly isolated fractional 
fluxes in a bulk material far from any boundaries requires 
a different explanation. 

This paper is organized as follows. Section HT1 describes 
the growth and preparation of the high-quality highly 



underdoped YBCO crystals and introduces the scanning 
Hall probe microscope. Section IIIII discusses Hall probe 
observations of partial vortices, and in Sec. II VI we model 
these partial vortices as kinked stacks of two-dimensional 
(2D) pancake vortices. Section [V] presents other exper- 
imentally observed properties, which are all consistent 
with the kinked vortex picture. Finally, in Sec. I VII we 
discuss the partial vortices in light of energy costs and 
pinning. The Appendix of this paper discusses the in- 
plane penetration depth extracted from fits to vortices 
in the YBCO. 



II. SAMPLES AND METHODS 

The YBa2Cu306+a; crystals, with x from 0.34 to 0.375, 
are grown with a self-flux method in BaZrOa crucibles de- 
tailed elsewhere. 8,9 After growth, the desired oxygen con- 
tent is set during a 900-930°C anneal in flowing oxygen, 
then oxygen inhomogeneities are removed during a 1-2 
week 570°C anneal in a small tube with YBCO ceramic 
at the same oxygen content. Initially, after quenching 
to 0°C, the crystals are non-superconducting, but an- 
nealing at room temperature allows the oxygen atoms 
to order into Ortho-II (every other chain empty) chain 
fragments whose increasing length provides the carrier 
doping in the Cu02 planes. 9 The superconducting tran- 
sition temperature (T c ) increases with room temperature 
annealing until saturation is reached after several weeks, 
giving final T c values of 5-20 K with bulk susceptibility 
transition widths (10%-90%) of less than 2 K£ This early 
generation of underdoped samples sometimes had a small 
(<2%) volume fraction of the 50-60 K T c phase of Ortho- 
II YBCO, as observed in magnetization measurements. 3,0 
During the room temperature annealing, a single crystal 
can be observed at a range of T c values. The platelet 
shaped crystals are about 1 mm x 1 mm wide, with their 
surface parallel to the a6-plane, and are typically 10- 
100 /im thick. As grown, the crystals have twinning 
boundaries, but they can be detwinned under uniaxial 
pressure at elevated temperatures. 



2 



Our most detailed observations of sub-<I>o par- 
tial vortices were made with a scanning Hall probe 
microscope , 11 i 12 ' 13 i 14 described in Ref. [ll|. Scanning Hall 
probe microscopy of single vortices is an established 
tcchniqu o 12 ' 13 i 14 first demonstrated by A. M. Chang 
et al.— Our Hall probe was made from GaAs/AlGaAs 
two-dimensional electron gas and had lithographic size 
0.5 /zm x 0.5 /im. The Hall probe measures the perpen- 
dicular magnetic field in the active area (with a constant 
offset). The Hall cross was covered by a thin film of 
gold which was grounded during operation. This gate 
prevented any stray electric charges on the sample sur- 
face from perturbing the Hall signal. For positioning in 
the z direction, the probe is mounted on the end of a 
thin aluminum diving board which forms a parallel plate 
capacitor with a copper pad underneath it. As the tip 
of the probe approaches the sample, we monitor the ca- 
pacitance and can ideally determine the location of the 
sample surface to within 10 nm. This touchdown pro- 
cedure is repeated at multiple locations within the scan 
area and then the probe is scanned in a plane just above 
the sample surface. The minimum height of the Hall 
cross active area above the sample surface is determined 
by the sample-probe alignment and for these measure- 
ments was 0.4 /im or larger due to geometric constraints. 
The lateral scan range of the microscope at 4 K is 60 /im 
and the sample can also be repositioned using xy stick- 
slip course motion. For improved signal-to-noise ratio, we 
averaged multiple images (having checked that the con- 
secutive images did not show changes). We sometimes 
used an 8 ^m x 8 ^m scanning superconducting quan- 
tum interference device (SQUID) with better flux sensi- 
tivity, but worse spatial resolution than the Hall probe. 
The cryostat was inside triple-layer mu-metal magnetic 
shielding with a residual field of less than 25 mG. 

The main results presented in this paper are from an 
8 ^m thick twinned YBa 2 Cu3 06.375 crystal imaged at 
eight stages during the room temperature oxygen order- 
ing annealing. The annealing took place in the micro- 
scope in a helium atmosphere. After 36 hours of anneal- 
ing, the crystal had T c ~ 5.1 K and transition width 
AT C ~ 3 K. Further annealing gave a range of T c values 
all having AT C < 1.5 K. The maximum measured T c was 
14.7 K (Fig. [1] inset). T c values were obtained in situ 
in an 8.3 mHz applied field of amplitude 0.20-0.25 Oe 
(Fig. Q]). The transitions are described as midpoint T c 's 
with full widths limited by the ~10% resolution of the 
susceptibility measurement. 

We also imaged flux in nine other similarly prepared 
YBa2Cu306+x crystals with x in the range 0.35-0.375 
and T c 's in the range 7-17 K. Each of these crystals 
was studied at only one T c value. Sub-<i>o flux features 
were seen in the three crystals with the lowest T C ~7K 
values lii The higher T c samples with T c > 11 K only 
showed flux consistent with conventional $o vortices. We 
also imaged vortices in Ortho-II (T c w 60 K) and near- 
optimally doped YBCO crystals and did not see any ev- 
idence of sub-$o partial vortices. 
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FIG. 1: Superconducting transitions in the YBa2Cu3C>6.375 
crystal measured by in situ magnetic susceptibility with the 
Hall probe in an applied 8.3 mHz field of amplitude 0.25 Oe 
for stages 1-2 and 0.20 Oe for stages 3-8. Stages offset by 
0.15 fi/G. Filled symbols are second measurements at dif- 
ferent locations. A Hall probe response of Rh = 0.115 ± 
0.015 fi/G indicates no measurable Meissner response, while 
fi/G indicates full shielding. Inset: Midpoint transition 
temperature (T c ) versus cumulative room temperature an- 
nealing time for the YBa2Cu3 06.375 crystal. Vertical bars 
indicate the full transition width, limited by our measurement 
resolution. 



III. FLUX IMAGES 

We saw over 100 sub-^o nvrx features in the 
YBa2Cu3 06.375 crystal while tuning T c from 5 to 15 K. 
We also observed more than 300 apparently full $o vor- 
tices, which were dominant for T c > 11.5 K. Typical 
Hall probe images are shown in Fig. [5] for a range of T c . 
These images were taken at low temperatures (T < T c /2) 
after field-cooling the crystal at about 3 K/min in a per- 
pendicular field. Though not necessary, the images were 
taken after turning off the applied field at low tempera- 
ture. Images taken before and after turning off the field 
looked identical. No flux features were observed when 
we cooled the sample in zero field. The images were not 
all taken at the same place on the crystal, because we 
occasionally used our course motion capability to move 
the sample in order to view nearby regions. 
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FIG. 3: Histogram of partial vortex peak field as a fraction of 
full vortex peak field in the YBa2Cu30e.375 crystal. Sample 
was cooled to low temperature in a perpendicular field before 
each image was acquired. Total bar heights show data from 
all anneal stages (all T c values). 



FIG. 2: (Color). Hall probe images of magnetic flux in the 
YBa2Cu3 06.375 crystal for increasing T c (stages 1-6). The 
probe's size was nominally 0.5 /im. The circular features with 
largest B z are likely straight pancake vortex stacks ( "full vor- 
tices"), while the dimmer and sometimes non-circular fea- 
tures are identified as segments of kinked pancake vortex 
stacks ("partial vortices"). The applied field H z while cooling 
through T c and the temperature T at which each image was 
acquired were (a) 0.15 Oe and 2.0 K, (b) 0.21 Oe and 2.3 K, 
(c) 0.21 Oe and 2.4 K, (d) 0.21 Oe and 2.2 K, (e) 0.43 Oe and 
4.3 K, and (f) 0.20 Oe and 4.1 K. A constant background 
has been subtracted from each image. For (f), low frequency 
telegraph noise in the Hall probe signal was subtracted from 
the raw images before averaging. 



The images in Fig. [2] show flux features that can be 
divided into two types. We identify the brightest fea- 
tures, which are close to circular, as conventional or "full" 
vortices. They carry total flux 4>o through the crys- 
tal surface, within experimental error. The full vortices 
increased in peak B z and decreased in width as T c in- 
creased, likely due to changes in the in-plane penetration 
depth (see Appendix). Other features have a smaller 
peak B z and appear either circular, elongated, or with 
tails. We call these features "partial vortices" because 
they carry less than <E>o of total flux through the surface. 
When the sample was cooled in a perpendicular field, par- 
tial vortices accounted for more than half of the observed 
flux features for T c < 11.5 K (stages 1-4), but dropped 
to less than 10% for T c > 14 K (stages 7 and 8). 11 We 
will show that the partial vortices can be explained by 
non-axial arrangements of 2D pancake vortices. 

The flux carried by the partial vortices was not re- 
stricted to discrete fractions of $o > as shown by a tally of 
the ratios of the peak B z of a partial vortex to the peak 
B z of a full vortex in the same image (Fig. [3]). The peak 
B z ratio for a partial vortex does not translate directly 
to the fraction of a $o through the surface, but nonethe- 



less, Fig. [3] indicates that the partial vortices occurred 
for a range of magnitudes. It is not feasible to directly 
tally the flux of each partial vortex without large errors 
because the field from nearby vortices in many cases in- 
terferes with the integration. 

A tempting hypothesis is that a partial vortex consists 
of a straight stack of pancake vortices, one in each layer, 
with some unusual mechanism by which the total flux 
is permitted to be less than <E>o- In this case, the field 
profile of a partial vortex above the sample surface could 
be calculated using the anisotropic London model (see 
Eq. (|A.1[) in the Appendix) except with a flux smaller 
than <I>o. If the in-plane penetration depth was assumed 
to be constant throughout the sample at each T c , then 
the peak B z ratio given in Fig. [3] would be equivalent 
to the fraction of a flux quantum carried by the partial 
vortex. However, we do not believe this hypothesis is the 
best explanation of our data. 

The main sources of error in the peak B z ratios were 
noise in the data and error in background determination, 
roughly 50 mG each. This gives an error in the peak B z 
ratios of ~5% at the higher T c values, and several times 
this for the lowest T c . Counts may also be missing at the 
ends of the histogram, since a partial vortex with close 
to full peak would likely be mistaken as a full vortex and 
those with very small peak field may have been lost in the 
noise. We omitted some flux features from the histogram 
if they could not be clearly identified as partial or full 
vortices. We also omitted whole images if they had high 
noise or did not contain a full vortex, the latter being 
particularly an issue for the lowest anneal stage where 
full vortices were rare. Only 21% of the images from 
stage 1 were included in the tally, while most images 
from subsequent stages were included. 
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IV. KINKED PANCAKE STACKS 

In this section we discuss a model in which the partial 
vortices result from kinked stacks of pancake vortices. We 
show that the model quantitatively describes the most 
circular partial vortices, and qualitatively describes the 
shapes and tails of the non-circular ones. 

In Ref. [3, Clem introduced the idea of 2D pancake vor- 
tices as the basic building blocks of 3D vortices in highly 
anisotropic superconductors consisting of weakly Joseph- 
son coupled layers. Even when the inter layer Joseph- 
son coupling is not negligible, the vortex structure can 
be described as a superposition of 2D pancake vortices 
and short sections of Josephson vortices (called "strings" ) 
connecting pancakes in adjacent layers^ Clemi also 
showed that for a vortex aligned along the c-axis, a 
straight stack of 2D pancake vortices gives the same re- 
sult as an ordinary 3D vortex in the anisotropic London 
model. However, Benkraouda and Clemi proposed that 
a tilted pancake stack may lower its energy by instead 
forming a kinked structure similar to the one shown in 
Fig. SJc), rather than maintaining a homogeneous tilt 
angle. Our observations of partial vortices suggest such 
configurations of kinks and short pancake stacks, which 
are stabilized by pinning effects. 

Compared to optimally doped YBCO, where 7 = 
Ac/A a b w 5.5±l^our highly underdoped Yl^CusOg+z 
(x ~ 0.35-0.375) crystals are much more anisotropic. 
Microwave measurements in similar samples found zero 
temperature values A c (0) pa 100 /im for T c pa 6 K, and 
A c (0) fa 40 fim for T c ps 15 K. 17 The zero tempera- 
ture in-plane penetration depth can be obtained from re- 
cent measurements^ of i? c i(0) in similar crystals which 
found that the power law H cl (0) = 0.366 T c 164 (Oe) fit 
H cl vs. T c data well for T c < 22 K. Using H cl = 
$o(ln(/t) + 0.5)/(47rA 2 b ) and k ps 40 (Ref. M) gives 
A Qb (0) = 1.00 /mi and 0.47 /mi for T c = 6 K and 15 K, 
respectively. These A a & values are in agreement with val- 
ues we obtained from fits to Hall probe vortex images 
from our variable T c YBa2Cu3C>6.375 crystal as discussed 
in the Appendix. Thus 7 ps 100 for T c pa 6 K and 7 pa 85 
for T c pa 15 K in these crystals. 

Since the 2D single layer screening length A = 
2A 2 lb /s pa 1 mm is greater than A c in our highly un- 
derdoped YBCO crystals (the bilayer spacing in YBCO 
is s = 1.17 nm), interlayer Josephson coupling is not 
negligible- and the notion of purely magnetically cou- 
pled pancakes is not entirely accurate. However, due 
to the large values of A c and 7, the additional attrac- 
tion between separated pancake stacks due to Josephson 
strings is smaller than the magnetic interaction, as will 
be discussed in Sec. IVI1 Thus the pancake vortex plus 
Josephson string picture should be at least qualitatively 
appropriate for these samples. 

To show that short pancake stacks can explain our ob- 
servations, we consider the magnetic field that a straight 
partial stack extending from Zi to Zf (zi < Zf < 0) gener- 
ates above the z — surface of a layered superconductor 



assumed to be much thicker than \ ab . It was shown in 
Ref. [TBI that the z component of the magnetic field at a 
height z above the surface and a radius r from the vortex 
axis is 

^Kb Jo QiQ + 1) 

(i) 

where Q = \J q 2 + A~ b 2 and the layer spacing s is much 
smaller than both z and A a (,. If the partial stack extends 
from \zi\ 3> \ a b to Zf = 0, Eq. fTJ gives the field of a 
conventional 3D vortex (Eq. (|A.ip ). Integrating Eq. ([1]) 
gives total flux 

* = $o(e*' /Aa * -e 2 '/ A ^) (2) 

through the z — surface of a superconducting half 
space for a partial vortex extending from z% to z/ . 15 i 20 
If an otherwise straight vortex stack has one kink at 
a depth z = —d, the total flux through the surface 
from the lower and upper partial stacks is <&$e~ d / > ' ab and 
<i> (l — e~ d / Aab ), respectively. The total flux is $0, as 
expected. 

To compare the kinked stack model with our observa- 
tions, we fit the two partial vortices labeled in Fig. 0Ja) 
to a model of a pancake stack with one kink. Figure |U(b) 
shows that partial vortex #1 has a narrower profile, so 
we chose it as the upper stack. Adding the field profiles 
(Eq. (Q])) of all the partial stacks in a kinked stack gives 
the field profile of a straight stack (full vortex). Figure 
HJb) shows that the cross section through a full vortex is 
similar to the sum of the cross sections through partial 
vortices #1 and #2. 

Using non-linear regression, we fit the portion of the 
image shown in Fig. H^d) inside the oval to a numerical 
approximation of Eq. (JTJ) with one expression each for the 
lower (—00 to —d) and upper (— d to 0) partial pancake 
stacks. Free fit parameters were the depth of the split, d, 
and the radial centers (xi,yi) and (£2,2/2) of each partial 
stack (Fig.[J)(c)). Fixed input values were z = 1.0 /jm and 
Kit = 0.65 /im. This X a b was obtained from a 2D fit to 
the full vortex labeled 'FV in Fig. 0|a) by the method 
described in the Appendix. Our fit gives a depth of the 
kink d — 0.36 /im and lateral displacement at the kink 
(kink length) p = y/fa - x 2 ) 2 + (yi - y 2 ) 2 = 5.3 /im. 
Figure |3Je) shows a 2D color plot of the fit and (f) shows 
cross-sections of the data and fit images along the dotted 
lines. Kinked or separated pancake vortex stacks is a 
plausible interpretation of our observations since, as Fig. 
@] shows, the model fits well to our data. 

This partial vortex pair was ideal for fitting because 
there appeared to be only a single kink within a few X a b 
of the surface and the partial stacks were well defined 
and circular. In most of our partial vortex images there 
were multiple kinks in a stack or several intermingled 
kinked stacks. In these cases, fitting the data would be 
more complicated. It was also common for partial vor- 
tices to have non-uniform shapes (see Fig. [^J) which do 
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FIG. 4: (Color). Fit to a kinked pancake vortex stack in the YBa2Cu30e.375 at anneal stage 6 (T c = 14.0 K) at T = 4 K. (a) 
Hall probe image containing partial vortices labeled #1 and #2. The 3 brightest features are full vortices, (b) Horizontal cross 
sections through the partial vortices, the sum of the partial vortices, and the full vortex labeled FV in (a). The solid line is 
from a 2D fit to the FV with z = 1.0 /jm as a fixed parameter and \ a i, — 0.65 /im as a fit parameter, (c) Sketch of the assumed 
geometry of the kinked pancake stack, not to scale, (d) Zoom in of image (a) with the color scale adapted to the signal range 
of the partial vortices, (e) Fit to the partial vortices within the white oval in (d) as discussed in the text, d — 0.36 /im and 
kink length is 5.3 /im. (f) Cross-sections through the data (circles) and fit (line) along the dashed lines in (d) and (e). 



not strictly agree with the model of one or a few kinks 
in a pancake stack. Elongated partial vortices and those 
with apparent tails, such as the partial vortex to the up- 
per left of #1 in Fig. 0Ja), could be the result of many 
closely spaced kinks, a tilt of a partial stack, or even a 
non-uniform staggering of pancakes from layer to layer. 
In each of these scenarios the displacement between ad- 
jacent pancakes would be smaller than the Hall probe 
spatial resolution. Though the kinked stack model can- 
not be used to fit irregularly shaped partial vortices with- 
out many free parameters, the underlying phenomenon is 
similar, with 2D pancake vortices playing a critical role. 



V. OTHER CHARACTERISTICS 

A number of observed properties of these partial vor- 
tices further substantiate the partial pancake stack in- 
terpretation as well as give further insight into their sta- 
bility and pinning. Partial vortices occurred in groups, 
preferred certain regions in the crystal, were more mobile 
than full vortices, and were more likely to be formed by 
cooling in a tilted field. 

The observed grouping of partial vortices is necessary 
for a kinked pancake stack, since all pancake vortices near 
the crystal surface carry flux through the surface which 



collectively adds to <&o. We observed partial vortices up 
to tens-of-microns away from others in a group, so even 
if a partial vortex appeared isolated, other segments may 
have been outside the image area. One caveat would be 
if there were subsurface vortex termination, which could 
occur for small sample size^i 

Partial vortices showed a tendency to prefer certain 
regions of the crystal, even after a room temperature 
annealing. This was especially noticeable for the later 
anneal stages (higher T c ) for which partial vortices were 
rare. This may indicate that kinks occurred preferen- 
tially in regions of the sample which were different from 
the bulk, perhaps with higher disorder, more pinning 
sites, or weaker superconductivity. At the lowest T c 
stage, we observed that partial vortices were more likely 
to pin where the tip of the Hall probe sat when the 
probe's active area was centered over the scan area. This 
was the location of the probe during cooldown through T c 
and also the place where the z approach was most often 
performed to determine the location of the sample. The 
vortices may have been preferentially attracted to the lo- 
cation of the probe's tip, or the repeated contact with 
the sample in that location could have created pinning 
sites. 

We also found that partial vortices were more mobile 
than apparent full vortices in the same samples. For 
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FIG. 5: (Color). Partial vortices coalesced during sample 
coarse motion (T = 2.4 K). (a) Hall probe image after cooling 
through T c = 8.6 K in H z = 0.21 Oe applied field. White 
circles identify two partial vortices, (b) After slight z coarse 
motion of the sample. Comparison of the images suggests 
that the partial vortices in (a) collapsed to one full vortex in 
(b). 



example, partial vortices sometimes moved or coalesced 
after stick-slip coarse motion of the sample holder in the 
xy or z directions. This was not observed for full vortices. 
Figured] shows Hall probe images before and after several 
ramps of the voltage on the z piezoelectric and slight 
z coarse motion. Two partial vortices of similar peak 
amplitude in Fig. [Ha) appeared to coalesce into one full 
vortex in Fig. GUb). Motion of partial vortices during 
coarse motion could be due to stray fields from the stick- 
slip high voltage pulses, which might create forces large 
enough to unpin some partial vortices. Once unpinned, a 
kinked stack could realign to a straight stack as favored 
by electromagnetic coupling of the pancakes and by any 
Josephson coupling. 

To encourage partial vortex formation, we cooled the 
sample in a magnetic field with a horizontal compo- 
nent to reduce the energy cost of a kink. With the 
YBa2Cu30e.375 crystal almost fully annealed with T c = 
14.4 K (stage 7), the sample was cooled through T c in 
an applied field H — H x x + H z z. The vortex arrange- 
ment did not change when the field was turned off at 4 K. 
The images in Fig. [5] show increased numbers of partial 
vortices with increased H x . The magnitude of H z deter- 
mined the density of flux observed in the images. The 
horizontal field may have caused pancakes within a stack 
to pin at large displacements with respect to each other 
as T was lowered through T c . Since the flux arrange- 
ment did not change when the field was turned off, the 
pinning must have been sufficiently strong to overcome 
the restoring forces favoring a straight stack. 

After cooling in a tilted field, we observed a change in 
the flux arrangement after the temperature was raised 
but still kept below T c . As shown in Fig. [7fa), the 
YBa2Cu3 06.375 sample with T c — 14.4 K (stage 7) was 
cooled to 3 K in a field H = 2.0x + 0.25z Oe, where H x 
was chosen to induce partial vortices. While at 3 K, the 
applied field was turned off and no change was observed. 
Then the sample was warmed to T = 6.6 K, cooled back 
to 3 K, and imaged again. This cycle was repeated several 
times with successively higher maximum T . After all cy- 




FIG. 6: (Color). Effect of an in-plane field on partial vortex 
formation in YBa 2 Cu 3 6 .375. T c = 14.4 K and T = 4 K 
for all the images. Field-cooled through T c in H z = 0.50 Oe 
and (a) H x = 0, (b) H x = 2.3 Oe. Field-cooled through 
T c in H z = 0.20 Oe and (c) H x = 0, (d) H x = 1.2 Oe, (e) 
H x = 2.3 Oe. Images acquired in zero applied field. The 
^-direction is horizontal and z is out of the page. 



cles, the flux arrangement changed. Vortices created by 
cooling through T c in H — 0.25z Oe did not show any mo- 
tion after similar thermal cycling (Fig. Efb)). These ob- 
servations indicate high metastability of the vortices for 
cooling in non-zero H Xl as there are many possible struc- 
tures consisting of pancake vortices and Josephson vor- 
tices in layered superconductors when cooled in a tilted 
field. 



VI. DISCUSSION 

The high-quality fit of the kinked stack model to our 
measured field profiles (Sec. IIV|) . as well as the other 
observed properties of the sub-<I>o features (Sec.|V]), sug- 
gest that we have observed kinked pancake vortex stacks 
in the highly underdoped YBa2Cu30g+a; crystals. The 
large anisotropy and large A c in these crystals suggest 
that it is reasonable to think of vortices as being com- 
posed of 2D pancake vortices. Strictly speaking, however, 
the Josephson coupling in these crystals is not negligi- 
ble compared to the magnetic coupling because A c is not 
much larger than A A The presence of Josephson coupling 
will lead to some distortion of the field profiles and to ad- 
ditional interaction energy, though a fully formed Joseph- 
son vortex will not exist along the kink since the kink 
length p is typically less than A c . Despite non- negligible 
Josephson coupling, the concept of pancake vortices is 
still qualitatively valid here and, as we have shown in 
Sec. IIV1 gives a good quantitative approximation of the 
field profiles measured with the Hall probe. 

Pinning is essential to explain our observations. We 
know that there are pinning sites in the YBCO crys- 
tals for all doping values, as is typical of type II super- 
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FIG. 7: (Color). Comparison of vortex movement after ele- 
vating the temperature for the YBa2Cua06.375 crystal cooled 
in tilted and perpendicular fields. T c — 14.4 K. (a) Field- 
cooled in H x = 2.0 Oe and H z = 0.25 Oe to T = 3 K. Image 
(i) was taken after the field was turned off. After thermal 
cycling sequentially to (ii) 6.6 K, (iii) 8.0 K, (iv) 9.8 K, and 
(v) 12.2 K. (b) The sample was field-cooled in H x = and 
H z = 0.25 Oe to T = 3 K. Image (i) was taken after the 
field was turned off. (ii) After thermal cycling to 11.8 K. 
T = 3.0 K for all images. 



conductors, because vortices remain in the sample af- 
ter the field is turned off below T c . In the absence of 
a field component in the a6-plane, magnetic coupling of 
the pancake vortices and any interlayer Joscphson cou- 
pling favors alignment of the pancakes along the c-axis. 
That we observe kinked pancake stacks in the highly un- 
derdoped YBCO even after the applied field is turned 
off indicates that pinning of the pancakes dominates the 
pancake vortex arrangement. It should be noted that the 
calculations in Sec. HVI are for kinks in otherwise straight 
pancake vortex stacks, but pinning may cause fluctua- 
tions in the pancake positions within a stack, distorting 
the field profiles^ 

To roughly quantify the strength of pinning required, 
we calculate the restoring force on the kinked structure. 
From Ref. 0, the energy required to deform a straight 
stack to a singly-kinked structure (in the limit of zero 
Josephson coupling and small layer spacing s) is 

E kink ( P ) = (^j J_[( e -P/A.»_i)+i n(p/Ao6 ) 

- Ei(-p/X ab ) + C], (3) 

where E^ is the exponential-integral function and C is 
Euler's constant. The restoring force on the kinked struc- 



ture with kink length p found from — dEkink/^P is 

(4) 

For a kinked stack like that in Fig. [3] with model param- 
eters p = 5.3 yitm and X a b = 0.65 /im, this restoring force 
is 80 fN. This result is for a kink deep in a crystal. 

We can also consider the additional cost of a kinked 
stack due to finite Josephson coupling. In the limit 
7s < p < A c , which is the case for all the kinked stacks we 
observed, Ref. [H gives the energy cost of a long Joseph- 
son string connecting the two partial stacks to be of 
the order Ej w (<j> /47r) 2 / 9(A a fcA c ) _1 . The corresponding 
restoring force is of the order Fj w — ($ /47r) 2 (A l bA c ) _1 . 
The A c and A a b values obtained from Refs. M and M 
give this force to be approximately 3 fN for T c = 6 K 
and 14 fN for T c — 15 K. Comparing this to .Fkink above, 
the added restoring force due to the Josephson string is 
somewhat smaller than that of the magnetic coupling. 
For a longer Josephson string, the energy cost of the 
string increases linearly with the length and the restoring 
force Fj remains constant. Thus, when the pinning force 
is large enough to compensate the restoring forces, there 
would be no theoretical limit to the maximum length of a 
metastable configuration. As a point of comparison with 
regards to pinning, measurements on a similar YBCO 
crystal with T c = 11 K have estimated the required force 
to unpin a full vortex to be ~0.5 pNi 2 ^ 

The calculated restoring forces would be much larger 
for a kinked stack in optimally doped YBCO. Using the 
values of X ab = 0.16 (Ref. S3) and 7 = 5.5 (Ref. [lj) 
for near-optimally doped YBCO, -Fkink has a maximum 
magnitude of 3 pN when p = 2X ab , or if p/X ab = 8.2 as it 
was for the highly underdoped calculations above, -Fkink 
would be 1 pN. In the limit of a long Josephson string, 
the Josephson coupling gives an additional restoring force 
of magnitude Fj = 2 pN. These forces are much larger 
than for our highly underdoped YBCO, which is not sur- 
prising since optimally doped YBCO has much smaller 
anisotropy and much smaller A c , thus we expect the vor- 
tex behavior at higher doping to be less two-dimensional. 
Indeed, we have never seen partial vortices in Ortho-II 
(T c « 60 K) or near-optimally doped YBCO crystals. 
Most of our measurements on these higher T c samples 
were after cooling in a perpendicular field. 

Pinning in these highly underdoped crystals must be 
strong enough to overcome the restoring force on the 
kinked pancake vortex stack. The origin of the pinning 
sites cannot be determined from our Hall probe images. 
They could be regions of oxygen inhomogeneity, twinning 
boundaries, lattice imperfections, or something else. Fu- 
ture studies of the pinning landscape in these crystals and 
in other cuprates is desirable. It is also not known how 
many of the pancake vortices must be pinned to support 
the kinked structure in the absence of an applied parallel 
field. However, the total force due to pinning must be of 
order 80 fN to compensate for the restoring force -Fkink- 
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The data show that the pinning becomes less sufficient 
for kinked vortex formation for the later anneal stages 
(higher T c ). As T c increases during the room tempera- 
ture annealing, X a b decreases and so the restoring force 
increases. Also as the crystal anneals the oxygen chain 
fragments get longer, and this may have an effect on the 
pinning. However, since the crystals have Ortho-II order- 
ing, there will still be many oxygen vacancies even when 
fully annealed. The number, spatial extent, and spatial 
distribution of vacancy clusters may all be changing at 
once during an anneal. An open question is whether oxy- 
gen inhomogeneities are necessary to see partial vortices 
in these samples. 

The behavior of the vortices when cooled in a magnetic 
field with a parallel component (Figs. [5] and [7]) also bears 
discussion. In a perpendicular field in the absence of pin- 
ning, a vortex should align parallel to the c-axis. When 
the field is tilted from the perpendicular, the ground state 
for a vortex depends on the parameters of the sample (see 
for example Ref. for vortex lattices). In our YBCO 
sample 7s < X a b and theory predicts a transition from a 
tilted to a parallel vortex lattice as the field approaches 
the afr-plane^ However, Benkraouda and Clem showed 
that beyond small angles a kinked structure can be en- 
ergetically preferable to an isolated tilted vortex, which 
itself is unstable beyond a tilt of 52°.— When our sample 
was cooled through T c in a tilted field we saw many more 
partial vortices, which we have suggested are kinked pan- 
cake vortex stacks. Since we did not see any changes in 
our images when the field was turned off at low temper- 
ature, pinning must have been sufficient to sustain the 
pancake vortex arrangement. However, when the sample 
was cycled to higher temperatures (still below T c ) the 
flux arrangement changed as shown in Fig. [7^a) . We hy- 
pothesize that at the higher temperatures the pancake 
vortices moved in an attempt to align along the c-axis, 
but in some cases became pinned again as kinked stacks 
upon cooling. In contrast, when the sample was cooled 
in a perpendicular field as in Fig.[7fb), the vortices were 
already pinned in their most stable state so no rearrange- 
ment occurred when cycled to higher temperature. 

The fact that the partial vortices formed even when the 
sample was cooled in a perpendicular applied field along 
the c-axis (within ~1°), especially for the lower T c values, 
further indicates the importance of pinning. Many of the 
sub-$o features are elongated or show tails, which could 
be due to staggering of adjacent pancakes or many unre- 
solvable kinks, indicating a complex pinning landscape. 
Recent improvements in the growth and preparation of 
these highly underdoped YBCO crystals have resulted in 
better homogeneity and thus perhaps less pinning com- 
pared to the earlier generation samples studied here. An 
ideal sample to study would be a YBa 2 Cu3 6 .33 3 crystal 
with Ortho-Ill inverse ordering (every third chain full). 
Such a highly underdoped crystal would have very few 
chain vacancies. Future work imaging flux in the next 
generation samples is desirable. 

A few other researchers have reported observations of 



kinked pancake vortex stacks. Grigorenko et al3 re- 
ported a Hall probe image of one "split" pancake vor- 
tex stack in a Bi 2 Sr2CaCu208+,5 (BSCCO) crystal with 
T c = 90 K. The kinked stack in that case was formed 
under a rapid change in magnetic field. Extensive vor- 
tex imaging has been done on BSCCO 2 ' 27 ' 28 and typi- 
cally kinked stacks such as ours have not been seen, in- 
stead combined or crossing lattice d 26 ' 29 of pancake vor- 
tex stacks and interlayer Josephson vortices are observed. 
Unlike in BSCCO, we would not expect a crossing lattice 
to appear as the applied field approaches the a6-plane 
in our highly underdoped YBCO crystals because the 
Josephson length 7s is smaller than \ a b^ 

Beleggia et al* 3 - observed dumbbell-like features con- 
sistent with kinked vortices in transmission electron mi- 
croscopy images of 300-400 nm thick films of optimally 
doped YBCO when the applied field was within 7° of par- 
allel to the film. Our work suggests that kinked vortices 
may form more readily at low doping and can form even 
in the absence of an applied parallel field. The partial 
vortices we observed after cooling in only a perpendicu- 
lar field could not have formed if sufficient pinning was 
not present. However, it is not known if pinning is al- 
ways required to form kinked stacks when a sample is in 
a continuously applied tilted field. 

In conclusion, we have observed sub-<I>o flux features 
in highly underdoped crystals of YBa2Cu306+a; with T c 
ranging from 5 to 15 K. These "partial vortices" are well 
described as segments of kinked stacks of 2D pancake 
vortices. The partial vortices were more mobile than un- 
kinked full vortices, formed more readily when cooled in 
a magnetic field with a horizontal component, and were 
seen most frequently for very low T c . Our observations 
provide a view of vortex behavior in the highly under- 
doped region of the YBCO phase diagram, showing that 
2D vortex behavior and an appropriate pinning landscape 
can produce complex flux features at the crystal surface 
that are distinct from conventional vortices. 
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APPENDIX: PENETRATION DEPTH 

We also used Hall probe images of vortices in the vari- 
able T c YBa2Cu3 O6.375 crystal to estimate the in-plane 
penetration depth A a t, which relates to the superfluid 
density n s oc A~ & 2 , as a function of T c and T. Our es- 
timates deviate from the well established linearity be- 
tween T c and n s (0) first suggested by Uemura et ai. 30 ' 31 
for higher doped cuprates. Our results supplement X a b 
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values obtained from H c i measurements in Ref. [l8j. The 
unknown height of the Hall probe above the sample gives 
large error bars on our X a b results, and certain caveats 
discussed below lead us to conservatively interpret our 
results as upper bounds on X a b- 

We chose a total of 40 different vortices in the 
YBa2Cu3 06.375 sample for fitting, all of which appeared 
to be well-isolated full vortices. Each vortex was fit with 
the anisotropic London model in the thick crystal limit 
with the a6-plane parallel to the surfacei 32 ' 33 i 34 ' 35 



B z (r, z) 



2ttA? 



dq 



qe qz J (gr) 
Q(Q + q) ' 



(A.l) 



where Q = y q 2 + X ab , r is the radial distance from the 
vortex axis, and z is the height above the sample surface. 
We integrated Eq. (|A.1[) at constant z over a 0.5 /jm di- 
ameter circular area representing the Hall probe. The 
results are insensitive to the exact probe size and shape. 
We fit the vortex images using non-linear regression to 
extract X a b with fixed z. Free parameters were the loca- 
tion of the vortex center, a constant offset in the magnetic 
field (due to the Hall probe), and X a b- The lengths z and 
X a b are strongly correlated and could not both be free 
parameters. 

In the scanning microscope z — zq + Az, where zq 
is the sample-probe distance when touching, and Az is 
controllable and for these measurements ranged from 0- 
0.16 /xm. Geometric constraints give a lower bound of 
zq > 0.4 /im. A very conservative upper bound of z$ < 
1.45 /im was obtained for this data set by fitting vortices 
at maximum T c with X a b = and z as a free parameter. 
A smaller upper bound of z$ < 1.3 /im was obtained 
by assuming X a b{0) is at least as large as in optimally 
doped YBCOi^i All vortices were fit with a range of zo 
values. Fit results are reported here with the typical 
value zq ~ 0.8 /im, with systematic error bars determined 
by fits with zo = 0.4 /xm and zq = 1.3 /im. 

For four vortices we also took Hall probe images while 
warming to investigate temperature dependence. The in- 
set of Fig. [5] shows X a b(T) at the minimum and maximum 
T c stages. Within our systematic and statistical errors, 
it is not possible to extract the details of X a b{T) at low 
temperatures. The penetration depth appears approxi- 
mately constant for temperatures below T c /2. Thus we 
approximate A af) (0) by A a6 (T min ) with T min ~2J? for 
T c < 12 K and T min - 4 K for higher T c . 

Our A a b(T m i n ) versus T c results are shown in Fig. [5] 
The dominant sources of error are the uncertainty in 
zq (which will be the same for all data points) and 
the calibration of the probe's Hall coefficient, Rh = 
0.115 ± 0.015 fi/G (which may fluctuate from cooldown 
to cooldown). Fits with the zq extremes show that the 
corresponding X a b error, indicated by the double arrow 
in the upper left of Fig. is roughly constant for all data 
points. The uncertainty in Rh translates to a X a b error of 
roughly ±8% of z + A a &. The in-plane penetration depth 
decreased as T c increased. For example, at T c — 6.5 K, 
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FIG. 8: Fit results for low-temperature X a b versus T c for 
the YBa2Cu30e.375 crystal, obtained from fits to scanning 
Hall probe images of 40 different vortices at eight T c stages 
of a room temperature annealing. The fitted X a b is an up- 
per bound on the true A a (,(0). Horizontal bars indicate the 
resolution-limited full width of the superconducting transition 
for each T c . Vertical error bars (shown only for the extreme 
data points) are from uncertainty in the probe calibration. 
A systematic error from uncertainty in the minimum sample- 
probe distance zo could shift the full data set by the extent 
indicated by the double arrow at T c = 4 K. Data points 
shown are with zo = 0.8 /im. For comparison, the dotted 
curve is from the fit to H c i(T c ) from Ref. [HI a s discussed in 
the text. Inset: Temperature dependence of the X a b results 
for two T c values. For each data set shown, X a b was obtained 
from fits to scanning Hall probe images of an individual vor- 
tex as T increased. The vortices disappeared before T c . The 
dashed lines indicate the maximum extent the data sets could 
be shifted due to systematic errors. 



the average value is X a b = 1.35 ± 0.18 1°'43 M m > while for 
the highest T c = 14.7 K it is A a6 = 0.57 ± 0.11 ±g;|g /im, 
where the quoted errors are from the uncertainties in Rh 
and zo, respectively. 

The in-plane penetration depth X a b is actually a com- 
bination of A a and A^, which are not equal in YBCO. The 
X a /Xb anisotropy has been measured to be slightly larger 
than unity in higher doped samples (see for example Ref. 
H). In-plane anisotropy would lead to a distortion of the 
otherwise circular field profile of a vortex aligned along 
the c-axis. In principle, our approach could be used to 
obtain separate values for A Q and A& if it was performed 
on a detwinned crystal. 

There are several caveats to our A a & measurements. 
First, it is possible that some of the full vortices identi- 
fied for fitting were actually partial vortices, especially at 
the lowest T c values where fewer full vortices were avail- 
able for comparison. If a partial vortex stack was fit with 
Eq. (jA.ip . the resulting X a b value would be falsely high. 
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At higher T c , many more full vortices of consistent ap- 
pearance were observed. However, it cannot be ruled out 
that the pancake vortices in the "straight" stacks were 
not completely axial but rather pinned in a staggered 
manner resulting in a more spread out vortex field pro- 
file. If this were the case, our fits would yield falsely high 
X a b values since Eq. (jA.ip is for a straight vortex. Thus 
the most conservative approach to our fit results is to 
take the values shown in Fig. [8] as upper bounds on the 
true values of A afc . 

For T c > 10 K, our X a b(T m i n ) values obtained from 
vortex fits are close to the values calculated from the low 
T c fit flci(O) = 0.366 T c 164 (Oe) from Ref. QJ, as shown 
by the dotted curve in Fig. [5) For T c < 10 K, our val- 
ues lie slightly above the curve. However, as discussed 



above, our values should conservatively be treated a s up - 
per bounds, so there is not a discrepancy with Ref. [l8l . 

Our Xat estimates, along with the results of Ref. [Lsl . 
give a larger T c for a given superfluid density in the 
highly underdoped regime than predicted by the Ue- 
mura relation T c oc n s (0). Other recent experiments have 
also shown deviations from this relation, such as Pereg- 
Barnea et al^ in higher doped YBCO crystals and Zuez 
et alM- in underdoped YBCO films. These data indicate 
that thermal phase fluctuations alone cannot explain the 
suppressed superfluid density in underdoped cuprates. 
Herbut and Case^ proposed that low temperature nodal 
quasiparticles and vortex fluctuations near T c can explain 
the observed nonlinearity between T c and n s (0). 
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